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Abstract 

An analytic equilibrium, the Toroidal Bessel Function Model, is used in conjunction with the 
gyrokinetic code GYRO to investigate the nature of microinstabilities in a reversed field pinch 
(RFP) plasma. The effect of the normalized electron plasma pressure (/?) on the characteristics of 
the microinstabilities is studied. A transition between an ion temperature gradient (ITG) driven 
mode and a microtearing mode as the dominant instability is found to occur at a (3 value of 
approximately 4.5%. Suppression of the ITG mode occurs as in the tokamak, through coupling 
to shear Alfven waves, with a critical j3 for stability higher than its tokamak equivalent due to 
a shorter parallel connection length. There is a steep dependence of the microtearing growth 
rate on temperature gradient suggesting high profile stiffness. There is evidence for a collisionless 
microtearing mode. The properties of this mode are investigated, and it is found that curvature 
drift plays an important role in the instability. 
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I. INTRODUCTION 



Confinement in the reversed field pinch (RFP) is dominated by global tearing modes under 
ordinary circumstances. Consequently, microscale drift-type fluctuations have received far 
less attention in the RFP than they have in the tokamak. However, operational modes such 
as pulsed poloidal current drive which applies an external force to flatten the current 
gradient, and the quasi-single-helicity state [2], whose quasi-spontaneous formation leads to a 
suppression of multiple helicity tearing modes, have gained success in reducing global tearing 
mode activity. In this regime, confinement time can improve significantly, and equilibrium 
gradients can steepen {3, 4]. This produces conditions in which microscale turbulence might 
emerge as a factor in confinement, just as it does in the tokamak. The experimental evidence 
for fluctuations that are independent of global tearing modes is still fragmentary but growing. 
The tearing modes couple to a broad cascade that delivers magnetic energy to small scales 
and simultaneously brings kinetic energy closer to equiparition as the wavenumber increases 
[ij], in keeping with MHD. However, at smaller scales, the kinetic energy eventually becomes 
greater than the magnetic energy. In this regime there are notable changes in the coherence 
and cross phase between magnetic field and density, likely associated with a different type 
of fluctuation. The observed radial structure indicates a standing wave pattern consistent 



with collisional shear Alfven waves or microtearing modes js]. There is also evidence that 
the temperature gradient in the edge is close to a critical gradient for temperature gradient 
driven microturbulence Q. 

After years of some uncertainty about what type of microinstability is most compatible 
with the high shear and low safety factor conditions of the RFP, numerical solutions of gy- 
rokinetic models in RFP geometry have shown that many modes familiar from the tokamak 
may also arise in the RFP js iOj. For example, ion temperature gradient (ITG) driven 
modes can dominate instability at zero (3, above its critical temperature gradient, albeit 
with a somewhat different character from their tokamak counterpart. ITG is a curvature 
drift resonant mode, and in the RFP it is the poloidal curvature that dominates. With 
bad curvature for all poloidal angles, the mode has a longer poloidal extent instead of a 
ballooning character [9[. The scale length of the critical temperature gradient for instability 
goes like the minor radius a instead of the major radius Rq. 

Analysis of instabilities that only arise for zero or very low values of (3 makes little sense 
for the RFP, where (3 values of 10% are not unusual. When increasing (3, it is reasonable to 



expect that the ITG mode is eventually stabilized - it is thus helpful to know at what fllJ^ 
stabilization occurs and how scales with various geometric and equilibrium parameters. 
It is also important to determine whether a new instability emerges at higher /3, and if so, 
at what critical (3 the instability arises and what characteristics apply to it. In the tokamak, 
the kinetic ballooning mode becomes unstable as (3 is increased beyond a certain threshold. 



In the RFP, it would not be surprising for a microtearing (MT) mode [Ill4l4j| to become un- 
stable for (3 values above the low-/3 ITG regime, provided the electron temperature gradient 
that drives the MT mode is comparable to the ion temperature gradient for ITG instability. 
The MT mode is a natural fluctuation for the RFP, if only because it is the small-scale 
extension of global tearing modes, which dominate confinement in ordinary discharges. The 
MT mode is not a current driven mode: current driven tearing modes require that a current 
gradient that varies on the scale of the minor radius a be larger than the flux discontinuity 
at the resistive layer, which goes like —2kg, ensuring that A' > 0, where A' is the standard 
parameter of the MHD tearing mode. With the binormal wavenumber kg large, microtearing 
modes have A' < 0. However, they can be driven by electron temperature gradients just 
like other drift modes, including trapped electron modes and electron temperature gradient 
driven modes. 

We show in this paper that under certain conditions, including high temperature gradients 
and weak density gradients, MT indeed supplants the ITG instability as the dominant mode 
above a critical value around 5% - 6%. In this transition the dominance of ITG at low 
(3 is assured by having an ion temperature gradient above the threshold a/Lxi ~ 3, where 
l/L Ts = VT S /T S is the temperature gradient scale length and s indicates the particle species. 
The dominance of microtearing at higher /3 is assured by having the electron temperature 
above a similar threshold [ajL^ e pa 3). We show that the coupling to Alfven waves that 
stabilizes the ITG mode as (3 increases does so at a higher (3 in the RFP than in the tokamak 
because of the shorter parallel connection length associated with the large poloidal field. The 
MT mode is identified on the basis of a mode structure that shows canonical tearing parity 
in the electromagnetic fields and a frequency in the electron direction. We examine scaling 
properties of the MT growth rate, including its threshold behavior, its scalings with the 
temperature ratio, and equilibrium quantities like the pinch parameter O. Growth rate 
scalings with respect to both r and G relate to magnetic shear scaling, which in the RFP 
is not an independent equilibrium parameter. This work follows a number of recent studies 
that suggest the microtearing mode may be important in both the tokamak and the RFP 



0, Il6l-ll9f . Here our analysis is specific to equilibria consistent with the Madison Symmetric 
Torus 2Gj. Moreover, in adapting gyrokinetic codes originally developed for the tokamak, 
we are careful to capture all of the effects arising from an equilibrium magnetic field with 
comparable poloidal and toroidal field components. This can be done using the toroidal 
Bessel function model in the GYRO code |21|, with RFP-appropriate representations of 



the curvature drift and parallel derivative 10] . We demonstrate here through benchmarking 



studies that quantitatively idential results can be obtained with corresponding modifications 
to the circular equilibrium model [22| of the Gene code 
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24]. 



An interesting aspect in recent gyrokinetic work relating to the MT mode is the obser- 



vation of instabi 
tokamak 
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ity in low-collisionality regimes. This is observed for simulations with both 
19] and RFP geometries js]. Theoretically, a series of papers specific to 



the tokamak in their approximations have collectively pointed to the conclusion that the 
microtearing mode is very sensitive to collisionality and should become stable for small 
collisionality 13KL51 |25| . One design of the current study has been to probe this mismatch 



between these theoretical predictions and the results of the aforementioned gyrokinetic simu- 
lations. A collisionality scan shows instability at low v, but behavior with other parameters 
suggests that there may be two branches of the instability, one for low collisionality and 
one for higher values more compatible with theoretical predictions. We specifically study 
the possibility, first suggested for RFP tearing modes in Ref. 26j, of an RFP MT branch 
that is enabled by the large electron curvature drift of the RFP. Artificial variation of the 
electron curvature drift strength shows that the growth rate diminishes toward zero when 
the curvature drift falls outside a certain range of values. 



This paper is organized as follows: First, the gyrokinetic model is introduced in Sec. [TT1 
Included in this section is a discussion of the geometric features of the RFP and how they 
are dealt with in GYRO, the gyrokinetic code used in the majority of these results. Then, in 
Sec. IIH] we show results from a (3 scan for a certain set of RFP parameters and discuss the 
instabilities seen at low and high f3. This is followed by a derivation of finite-/? suppression 
of ITG in Sec. II II Al Next, Sec. IIII Bl contains a detailed characterization of the RFP mi- 
crotearing mode, including an investigation of growth rate dependencies on parameters such 
as collisionality, temperature gradient, and curvature. All of these results are summarized 
in Sec. ED 
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II. GYROKINETIC MODELING 



A comprehensive method by which to investigate microinstabilities in toroidal plasma 
devices is the gyrokinetic framework, in which the fast gyromotion has been removed from 
the kinetic equations through a gyrophase averaging procedure [271 ]. reducing the problem 
from the original 6D phase space to a more computationally tractable 5D. 



In this work we present results from the initial- value gyrokinetic code GYRO [2l|. This 
code provides a numerical solution of the gyrokinetic Vlasov-Maxwell system of equations. 
We here restrict ourselves to linear simulations in a radially local (flux-tube) domain, in 
which the main quantities of interest are the growth rates and frequencies of the most un- 
stable mode at a particular wavenumber kop s , where p s = {T e /mi) x ' 2 mi/eB is the ion sound 
gyroradius. This linear analysis provides an initial characterization of microinstabilities in 
the Madison Symmetric Torus (MST) |28|]. An expanded nonlinear analysis is currently in 
progress, and this nonlinear work, in conjunction with measurements from up and coming 
diagnostics on MST, will provide a unique opportunity for validation of the gyrokinetic 
results. 

In order to properly capture the effects of the magnetic geometry of the RFP, a number 
of modifications have been made to the code. These changes concern both the magnetic 
equilibrium used and the expression of certain geometric operators. A more detailed expla- 



nation of the changes and their effect on linear growth rates is given in Ref. lOj, and only 
an outline will be given here. 

The fundamental change that has been made is the inclusion of the Toroidal Bessel 
Function Model (TBFM) as the equilibrium model used in GYRO. The standard equilibrium 
model in the RFP is the Bessel Function Model (BFM) \B e = B Ji(Ar), = B J (Xr)] 
which is the solution of the force free condition V x B = AB in cylindrical geometry. The 
TBFM is an extension of the BFM to toroidal geometry, and since it is derived from the 
Grad-Shafranov equation, it is able to incorporate finite pressure effects. Flux surfaces are 

assumed to be circular, and to lowest order in (3 the magnetic fields are given as: 

Bp^er/a) B J (2Or/a) 

1 + r cos r cos 9/Rq 

where O = (Bg) wa11 / (B^)™ 1 is the RFP pinch parameter. In these simulations, O is an 
important input parameter, from which the safety factor go and shear s are calculated self- 
consistently. The magnetic geometry in the local simulations is thus completely determined 
by O and the normalized radius r /a of the simulation domain. 



In addition to the background magnetic equilibrium, certain of the gyrokinetic operators 
are generalized to the RFP geometry. These generalizations stem from characteristics of the 
RFP equilibrium that differentiate it from the tokamak. Primarily, the poloidal magnetic 
field in the RFP is much stronger relative to the toroidal field than in the tokamak, and 
in the case where the normalized radius of the simulation domain is r /a ~ 0.5 (a typical 
region of interest in simulations), they are roughly the same order of magnitude. Under 
these circumstances, the common tokamak approximation B ~ cannot be made, and the 
more general form B = B^l + (e t / qo) 2 ) 1 / 2 must be used, where B^ is the toroidal magnetic 
field and e t = r/R is the inverse aspect ratio. This difference affects operators such as the 
curvature drift frequency and parallel transit operators, the latter of which in its general 
form will pick up an additional factor, 



Besides being incorporated into the code, these geometric modifications must also be taken 
into account in analytic estimates whenever generalizing a tokamak analysis to an RFP envi- 
ronment. The next section provides an example of this in the context of finite-/? suppression 



The results presented below were performed using a set of parameters chosen to represent 
the experimental conditions in MST. This device has a major radius of Rq = 1.5 m and 
a minor radius of a = 0.5 m, yielding an inverse aspect ratio of e t = 1/3. We look at 
the radial location r^/a = 0.5, with, unless otherwise stated, the other parameters being: 
go = 0.186,5 = -0.716,9 = 1.35, a/L n = 0.58, a/Lx = 5.0, Ti/T e = 0A,u(a/c s ) = 0. The 
collisional frequency u, in particular, plays an important role in the dynamics of the MT 
mode, which will be discussed in Sec. IIIIBI 

Now we look at the effect of (3 in determining the nature of the dominant linear instability. 

III. BETA SCAN 

A scan over /3 from to 10% was performed to determine the variety of modes that might 
be dominant across this range. The results can be seen in Fig. [TJ It should be noted that for 
this scan, the pinch parameter is kept constant and does not vary self-consistently with /3. 
At low 0, the dominant instability is identified as an electrostatic ITG mode. As (3 increases 
this mode is suppressed and eventually overtaken by a MT mode at a (3 of approximately 



b- V ~ k 



1 



1 



(2) 




of ITG. 
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4.5%. The transition can be seen most clearly in the real frequency plot. Here, the frequency 
of the ITG mode is in the ion direction (negative sign in this convention) and that of the 
MT mode is in the electron direction (positive sign). 

These results have been confirmed with the circular equilibrium in the Gene code for 
the (3 = ITG regime (see Fig. [2]). Since gyrokinetic simulations are a relatively new topic 
of study for the RFP, it is important to compare results from separate codes with slightly 
different handling of geometric quantities. In order to attain agreement between the codes 
it was necessary to generalize Gene's circular equilibrium by including two additional terms 
to d-B/dr: a term proportional to s, and a term accounting for the r dependence of B^. 

One interesting feature of the low-/? versus high-/? instabilities is the range of scales at 
which these modes are unstable. This difference can be seen in Fig. |3j At low /?, the ITG 
mode ranges from kgp s = 0.1 to kgp s = 0.9, achieving a peak near 0.5. Initially, the MT 
mode arises at these same scales, with a peak at roughly the same value, though with a 
slightly broader range. As (3 increases, however, the peak of the microtearing mode shifts 
to higher k$. This can be seen in Fig. [31 where it is observed that at a /3 = 9% the MT 
mode not only reaches higher growth rates, but does so over a much larger range of scales, 
peaking at a value of k e p s w 1.5. 

The range of /? simulated covers different modes of operation of MST. A standard dis- 
charge will have values of (3 ~ 4 — 5%, a range that means ITG and microtearing may both 
be equally strong. Improved confinement pulsed poloidal current drive (PPCD) discharges, 
on the other hand, may achieve (3 values of 9% or higher 29[, in which case microtearing 
may be the dominant mode. Runs with [3 = 9% will be the focus of Sec. IIIIBI In the next 
section, we discuss the observed suppression of ITG with increasing (3. 



A. ITG beta suppression 



Finite-/? suppression of ITG (a linear effect which may be amplified nonlinear ly 30]) has 



been a topic of study in the context of tokamaks [3l| , and that analysis will be applied here 
to ITG in the RFP. In the RFP, ITG growth rates can still be quite strong at values of /? 
where tokamak ITG is typically stable 21], as can be seen in Fig. [TJ In fact, ITG remains 
unstable past (3 = 5% (though subdominant to microtearing) and may not stabilize until 
/? ~ 10%. 



Our analysis of the finite-/? suppression of the ITG mode follows that of Hirose [31], 




FIG. 1. Growth rate and mode frequency plotted as a function of (3 for kgp s = 0.372. A transition 
of the dominant mode from ITG to MT occurs at (3 ~ 5% 




FIG. 2. (Color online) A comparison of Gene (red +'s) and GYRO (green x's) results for /3 = 
in the RFP. Other parameters are given in the text. The two codes show good agreement. The 
difference in real frequency at high kg may be a results of incomplete convergence. 
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FIG. 3. Growth rate and frequency as a function of wavenumber kgp s for two different values of 
j3. ITG is seen to be dominant at (3 e = 1%, microtearing at j3 e = 9% 
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FIG. 4. (Color online) Growth rate plotted against temperature gradient for k$p s = 0.372 in the 
case of MT (red squares) and ITG (green circles). Both instabilities have a threshold around 
q/Lt ~ 3.5 — 4.0, for their respective driving gradients. 
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making the appropriate modifications for the RFP geometry (see Sec. [TTJ) . As discussed, 
these modifications are due to the different strengths and scale lengths of the magnetic 
field. In the tokamak, the scale length of magnetic field variation is proportional to the 
major radius, 1/Lb — VB/B ~ 1/Rq, while in the RFP, the appropriate scale is the minor 
radius, VB/B ~ 1/a. As was mentioned above, the parallel derivative term also needs 
to be modified to account for the equivalent strengths of the poloidal and toroidal fields. 
Therefore, the modified terms will take the forms fcn = l/(qoRo(l + (et/qo) 2 ) 1 ^ 2 ) and the 
curvature drift uoe = 2cT e (VB x B) • k/e£> 3 ~ I/Lb, where c is the speed of light and e is 
the fundamental charge. 

The following discussion is based on a fluid analysis of the ITG instability. Before in- 
cluding the finite-/? effects, the ion and electron densities are found to be, respectively, 



(u + 5u Di /3) 2 - 10u4/9 



n i = 777TT-. /o \ 9 m„2 /n » ( 3 ) 



n e = —n , (4) 

where u* e = cT e (Vlnno x B) ■ ~k/eB 2 ~ 1/L n and rji = dhxTi/dhxnQ. 

The expression for the electron density is modified by the consideration of finite-/? effects. 
This is done by taking into account perpendicular magnetic field perturbations, or, equiv- 
alently, perturbations to the parallel magnetic vector potential: B^ = V x Ay. Including 
such perturbations in the parallel momentum balance of electrons will result in the electron 
density taking the form 



Using Ampere's law and the quasineutrality condition, we arrive at the following relation, 



P 



k 2 L n Lg 



2% + - (1 + 2e n ) + 1 + rj e 



T 2 T 



ck\\ k 2 , V T 



1 + 2— $, (6) 



where r = T e /Tj, kn e = (47mne 2 /T e ) 1 / 2 , r] e = d In T e /d In n , k± = k • Bj_ ps ke, and 



e n = L n /Ls- Then, as in Ref. 3l|, we are able to derive the stability condition, with certain 



terms adjusted to account for the RFP generalizations, 



s e 2 r 2 



P ~ (l + (e t /q y)q 2 [(T + 2E n )(T + l) + T 2 r ]e y (?) 
Choosing parameters similar to those used for the growth rates plotted in Fig. [TJ (r = 
2.5, e n = 1/. 58, Tj e = 8.6, go = 1/5, et = 1/3), the above expression yields a critical (3 of 
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approximately 10.6%. This agrees well with an estimate of fllJ^ from Fig. [T] based on 
extrapolation of the ITG growth rate for > 5%. 

The criterion given above is similar to that of Ref . {jy] , but for RFP parameters it yields 
a higher critical /3 than is seen in tokamaks. This is due primarily to the smaller go and 
shorter connection length that is a result of the equivalent strengths of the poloidal and 
toroidal fields in the RFP. 

We now turn our attention to the instability observed at higher /3, the microtearing mode. 



B. Microtearing Modes 

The dominant instability at higher values of (3 is identified as a MT mode. These modes 
are characterized by having tearing parity in the parallel direction. This appears as odd 
parity in the electrostatic potential $ and even parity in the magnetic potential An, as seen 
in the eigenmode structure in Fig. [5j This structure, including the small amplitude features 
occurring every 2n (and attributed to poloidal variation) is similar to that seen in other 
devices, including MAST 1(| and RFX-mod [8j. Several parameters scans were performed 
in order to better characterize the observed MT mode. These scans were performed at a /3 
value of 9%, chosen to lie in the range for PPCD discharges in MST. 

The MT mode is known to be driven by the electron temperature gradient (rather than the 
current gradient, as is the case for the large scale tearing mode). Thus, a strong dependence 
of the growth rate on Ly e is expected, and this is seen in Fig. IIII Al Here the mode can be 
seen to require the threshold gradient of approximately a/L^e ~ 3.5, which is higher than 
that reported in 8j ((a/Lr e ) crit ~ 2). This threshold falls at a similar value as is observed for 
the ion temperature gradient threshold for the ITG mode. However, it can be seen that the 
MT growth rate rises much more steeply with temperature. Such a strong dependence can 
be expected to lead to profile stiffness, fixing the experimental gradient near the threshold 
for instability. Additionally, nonlinear simulations in the spherical tokamak have revealed a 
nonlinear upshift in the effective gradient threshold as compared to linear simulations jlo| . 
and the same effect can be expected to occur in the RFP. 

In Fig. [6] can be seen the effect of varying the temperature ratio Tj/T e . MT shows a 
strong increase of the growth rate with this ratio. There is an opposite dependence in ITG 
(see plot b), where the growth rate decreases as Tj/T e increases. This is consistent with 
MT's electron temperature gradient drive. 
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FIG. 5. (Color online) Eigenmode structure for the MT mode in electrostatic potential $ and 
magnetic potential An with both real (green dashed curve) and imaginary (red solid curve) compo- 
nents. The fields are plotted against the magnetic- field following ballooning angle 8*. This mode 
displays tearing parity, which is recognized as even parity in An and odd parity in <I>. 
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FIG. 6. (Color online) Growth rate plotted against the temperature ratio Ti/T e for kgp s = 1.488. 
Shown are MT (a) and ITG (b). 
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Since collisions are expected to play an important role in the mechanism for instability 
of MT 12|, ll4| , a scan was performed over the collision frequency, the results of which are 
given in Fig. [3 The collision operator used in GYRO is a pitch-angle scattering operator 



c(/) = ^| (1 _ a |, ( 8) 

where £ is the pitch angle and ^ e (e) = v[Z c r + H(e 1 ^ 2 )]/e 3 ^ 2 , with Z cS the effective nuclear 
charge, H(x) = exp[— x 2 ]/x-k 1 / 2 + (1 — l/2x 2 )erf(x) and e = E e /T e = m e v 2 /2T e . In these 
expressions, v is the control parameter in the simulations. 

A notable feature of Fig. [7] is the appearance of what seem to be two separate regimes 
of MT, with a transition between the two occuring roughly around v ~ 0.1 — 1. At v ~ 1 
the growth rate achieves a peak and then falls off for higher collisionality. At lower values 
of v the growth rate flattens, remaining finite in the limit of zero collisionality (although 
GYRO uses an upwind differencing scheme, which may introduce collisional effects). It 
should also be noted that the real frequency scales linearly with v above v « 1. Previous 



analytic work in tokamaks (see Refs. 12|, ll4() has concluded that MT should be stable in 



the collisionless limit, though some evidence contradicting this conclusion has been seen in 



gyrokinetic simulations in tokamak |T7j and spherical tokamak [16j geometries. 

The collisional dependence of this mode was investigated at several different radii (varying 
only q and s in correspondence with the TBFM and keeping all other parameters fixed). 
The results of this can be seen in Fig. [HJ Importantly, the wavenumber spectrum (a) shows 
a general stabilization of the mode at larger radius, which might be attributed to the larger 
shear at these locations. The collisionality scan (b) also shows interesting behavior: in 
particular, the growth rate at low v is much more affected by increased radius than at 
v ~ 1, and at r/a = 0.6 the mode is completely stabilized in the collisionless limit. This 
behaviour may suggest that there are two distinct varieties of MT modes in these simulations 
- one at low collisionality and one at moderate collisionality - that have somewhat different 
physics behind their drive mechanisms. The effects of increasing the pinch parameter O are 
similar to the effects of increasing the radius. This is to be expected from the TBFM, due 
to the same parametric dependence on r and O in that model. 

The physics of the drive mechanism was investigated further by looking at the role of 
curvature drift in the instability. In this study, a scalar factor a was placed in front of 
the curvature drift term in GYRO and varied to change its relative strength; a value of 
a = 1 corresponds to the physical curvature drift. The results of this are seen in Fig. [9] 
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for v = 0.001 (a) and v = 1.0 (b). In both cases, the instability is strongest for a ~ 1 
and falls off for values much lower or higher than this. This behavior is consistent with 
the magnetic curvature drift instability derived by Finn and Drake (see Ref. 26|), in which 
they describe a semi-collisional drift-tearing mode unstable in the presence of curvature 
drift and an electron temperature gradient. Ref. 26j is formulated in a cylindrical RFP 
equilibrium, and they note that this instability is strong only when wo ~ a;*, a characteristic 
that makes it more relevant to the RFP, with its associated stronger curvature drifts, than 
in a standard tokamak. Spherical tokamaks might also be expected to be susceptible to the 



curvature drift instability. We also note that although Ref. [26] uses fluid theory in the semi- 
collisional regime, collisions do not play an explicit role in the instability, and it is plausible 
that this mode may arise in the collisionless limit with proper inclusion of curvature drifts. 



IV. CONCLUSION 

The linear characteristics of microinstabilities in RFP geometry were investigated using 
the gyrokinetic code GYRO. It was found that ITG is the dominant gyroscale instability 
at low (3 for the parameters used in this work, with an increase in /3 causing a transition 
to MT. An analysis of the finite-/? suppression of ITG demonstrated both qualitative and 
quanititative agreement with the simulation results. This analysis suggests that the sup- 
pression mechanism in the RFP is the same as in the tokamak, with a larger critical that 
arises predominantly from a shorter parallel connection length. 

The MT mode was further investigated by performing a variety of parameter scans. The 
characteristics of the mode were in agreement with previous simulation work in other devices. 
Evidence for a collisionless MT mode was seen, though further work remains to be done as 
to be determine the nature of the physics behind this mode. There is some evidence that 
curvature drift plays an important role in the instability, and it may be described by the 
magnetic curvature drift instability of Finn and Drake 26]. 

The results of this work regarding implications for the Madison Symmetric Torus suggest 
that MT may indeed be playing a role in the gyroscale dynamics. Standard discharges lie in 
the range of for which both ITG and MT may be excited, while PPCD discharges achieve 
(3 values high enough so that ITG may be sufficiently suppressed and only MT is important. 
For these conditions and a large enough temperature gradient, the strength of MT and 
range of scales at which it exists suggest that it should be making important contributions 
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FIG. 7. MT growth rate and frequency plotted against the collisional frequency v for kgp s = 
1.488. There appear to be two distinct regmes: a region of constant growth rate and constant real 
frequency at low v and a separate region at v with a peak in growth rate and a real frequency that 
scales linearly with v. 
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FIG. 8. (Color online) MT wavenumber spectrum (a) and collisionality scan (b) for different values 
of r/a. The corresponding values of qo and shear are - r/a = 0.4: qo = 0.209, s = —0.382 (red solid 
curve); r/a = 0.5: qo = 0.186, s = —0.716 (green dashed curve); r/a = 0.6: qo = 0.155, s = —1.344 
(blue dotted curve). There is stabilization of MT with increasing radius r/a, especially prevalent 
at low v. Increased radius coincides with increased shear, which may play a role in stabilization. 
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FIG. 9. (Color online) The role of curvature drift in the MT instability. The parameter a is a 
factor regulating the strength of the curvature drift in the code. Shown are v = 0.001 (a) and 
v = 1.0 (b), as well as kep s = 1.488 (red solid curve) and kgp s = 0.186. This behavior is similar to 
that seen in Ref. [26!]. The points at low a represent a separate mode that has not been studied in 
detail. 
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to gyroscale turbulence. Let us note, however, that for alternate sets of parameters other 
modes, such as the trapped electron mode or kinetic ballooning mode, may emerge and 
play a role in gyroscale dynamics. Further work must be completed to determine the critical 
parameters required for these various instabilities and the effect of the background equilibria 
on their growth rates and transport fluxes. This task, as well as measuring and characterizing 
these modes experimentally, is an important area of future work. 
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